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We report a detailed study of Eulerian and Lagrangian statistics from high resolution Di- 
rect Numerical Simulations of isotropic weakly compressible turbulence. Reynolds num- 
ber at the Taylor microscale is estimated to be around 600. Eulerian and Lagrangian 
statistics is evaluated over a huge data set, made by ISSG'^ spatial collocation points and 
by 16 million particles, followed for about one large-scale eddy turn over time. We present 
data for Eulerian and Lagrangian Structure functions up to ten order. We analyse the lo- 
cal scaling properties in the inertial range and in the viscous range. Eulerian results show 
a good superposition with previous data. Lagrangian statistics is different from existing 
experimental and numerical results, for moments of sixth order and higher. We interpret 
this in terms of a possible contamination from viscous scale affecting the estimate of the 
scaling properties in previous studies. We show that a simple bridge relation based on 
Multifractal theory is able to connect scaling properties of both Eulerian and Lagrangian 
observables, provided that the small differences between intermittency of transverse and 
longitudinal Eulerian structure functions are properly considered. 



1. Introduction 

In the last few years many interesting and remarkabe results have been obtained by in- 
vestigating the statistic al properties of Lagrangian pa rticles advected by a turbulent flow 
-for a recent review see Toschi fc Bodenschat 3 ()2009l) . For Lagrangian particles we mean 
simple point like ideal tracers, whose instantaneous velocity coincides with the local Eule- 
rian velocity field, u(x, t): X(i) = u(X(i), t). Beside being relevant in many applications, 
where transport and/or aggregation of particles is important, the study of fully developed 
turbulence in the Lagrangian framework has opened important new directions of scientific 
investigations. In particular, small scale intermittency of turbulent flows and dissipation 
range can be probed m ore efficiently in the Lagrangian framework as shown experimen- 
tally s ii ice the works of iLa Porta etall (120011 ) ; 'Mord ant et"aI1 (l200ll.l2006h:IOtt fc MannI 
( 2000l ): lBerg et al.l (|2006h : iBourgoin et all (|2006 ): .Luthi et all (|2005h : iBerg et all ^200% 
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and numerically in 



Biferale et all (l2005l): lY cung ct al.' feO Od): iHomann et all (|2007f ): 



Mordant et"aI1 (|2006l ). (see also lBiferale et ali (,2008,) and iArneodo et alJ (|2008l ) for com- 



parison between experimental and numerical Lagrangian data). Moreover, Lagrangian 
and Eulerian measurement (those on the reference frame of the labora tory) sh o uld b e 
intimately sta t istica lly connected as pointed out for the first time in Borgad ( 2003 ) 



Boffetta et al. 1 (l2002h .' the hope is to learn about the former from the study of the latter 



and viceversa. 

In this paper we present some new results concerning the statistical properties of the 
Lagrangian and Eulerian velocity field in weakly compressible turbulent fiow. The flow 
is the outcome of a high resolution numerical simulations described in the next section. 
Our motivation is threefold: due to the high spatial resolution and the huge particle 
statistics (16 millions of particles) we are able to extract the statistics of the Lagrangian 
velocity field with high accuracy and relatively small error bars. From this point of view 
we arc able to assess, for the first time at this Reynolds numbers, statistical fluctuations 
in the Lagrangian domain as intense as those given by moments up to order 10. The 
second purpose is to understand whether the effect of the weak compressibility changes 
the statistical properties of Eulerian and Lagrangian velocity field with respect to pre- 
vious numerical investigations of incompressible isotropic t urbulent flow s . Pre liminary 
results concerning this point have been already published in iB enzi et al. I (|2008l) : here a 



more detailed analysis on new observables is reported. In particular, we are interested in 
understanding the effect of intermittency in the inertial range for Lagrangian particles. 
Some results concerning dissipative scales will be also presented. Finally, our third moti- 
vation concerns with the statistical link between Lagrangian and Eulerian measurements. 
Concerning this point we show that a simple model, based on the multifractal theory, 
is able to translate between the two ensembles with a quantitative agreemen t extend- 
ing ov er three decades from inside dissipative scales to the integral scale (jArncodo et al.l 



Our main tool is the numerical computation of velocity increments over a time lag t, 
the so-called Lagrangian structure functions: 

Sl'Hr) = {[v,it + r)-v,{m = {{SrvOn, (l-l) 

where i = x,y,z are the three velocity components along a particle trajectory, Vi{t) = 
Ui{x{t),t) and the average is defined over the ensemble of particle trajectories. As sta- 
tionary and homogeneity is assumed, moments of velocity increments only depend on 
the time lag t. When isotropy is also valid, all components must be symmetric and we 
will drop the dependency on the spatial index. In the inertial range, for time lags smaller 
than the integral time and larger than the Kolmogorov time, r,, <^ r <C Tl, non- linear 
energy transfer governs the dynamics. From a dimensional viewpoint, only the scale t and 
the energy transfer e should enter the structure functions. The only admissible choice in 
isotropic statistics is S^p^ (t) {ery/^ , but it d oes not take into account the fluctuating 
nature of dissipation. Many empirical studies ( Toschi fc Bodenschatj ( 2009( )) have in- 



deed shown that the tail of the probability density functions of StV become increasingly 
non-Gaussian at decreasing t/Tl, leading to intermittency and to anomalous scaling 
exponents, meaning a breakdown of the dimensional law, i.e. 

S(p\t)^t^''' (1.2) 

with ^(P) ^ p/2. 

As we will review and test in details later, the Lagrangian measurements can be linked 
to Eulerian Structure Function (ESF). In isotropic statistics we may have two differ- 
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ent Eulerian increments, longitudinal and transverse (see iFrischI ( 19951 ) for a textbook 
introduction to scaling in homogeneous and isotropic turbulence): 



4P)(r)^([u(x + r)-u(x)].fn 
4P)(r) = (|u(x + rT)-u(x)|f), 



• u = 0. 



(1.3) 



It is well established, ([Arneodo et al. I (|l996l) : iFrischI (|1995D ) that also Eulerian statistics 
have anomalous scaling, for <C r <C L: 



5f (r) 
S^^\r) 



(1.4) 



From a theoretical point of view, in isotropic turbulence one would expect that longitudi- 



Ap) 



see 



Biferale fc Procaccia 



nal an d transverse fluctuations have the same scaling, (j^ ^ 

(|2005l) : or a review on anisotropic t urbulence). This is not what was observed in many ex- 

perini ei ital and numerical a nalysis (Boratav fc Pelz ( 19971): van der Wateij (ll999l ): IShen fc WarhaftI 

~a)) 



(I2OO2I ): IChen et all (Il997l ): IZhou fc Antonia 



can be due to finite Re ynolds numbers (iHill 



( 20001): iDhruva et al 



( 20011): iHe et all (|1998I) ) or to some remnant 



( 19971 )'). Discrepancy 



small-scale anisotropy (jBiferale fc Procaccia Still, even large Reynolds number 



DNS wit h isotropic forc i ng sh o w some persi s tent d ifferences between (^'^ and Ct'-' as 



shown in Ishihara et al. I (|2009l) : iGotoh et all (|2002h : the same discrepancy is quantita- 
tively confirmed in the analysis here reported. Whether this effect will persist at higher 
Reynolds number remains an important open question for further experimental and nu- 
merical investigations. 

Recently, it has been argued that there must exists a relation between L agrangian and 



Euleri a n exponents, on the basis of a common multifractal description (jBoffetta et al 



(|2002l ): lBiferale et all (120051) ). This relation between Lagrangian and Eulerian statistical 
properties is highly non trivial and it is worth to understand how accurately can be 
verified in the existing DNS and laboratory experiments. It has b een proved t o be very 
efficient to predict the pdf of Lagrangian a cceleration (Bif e rale et al. ( 2004 )) and of 
Lagrangian Structure Functions of low order ( Arneodo et al.l ( 20081 )). 

The paper is organized as follows. In section ^ we discuss the numerical simulations 
and the data set. In section ([3]) we discuss the computation of the Eulerian structure 
functions. In section ([4]) we report the analysis of Lagrangian structure functions and we 
discuss the physical results obtained by our analysis. In section ([5]) we provide a theoret- 
ical interpretation of our findings by employing the multifractal framework. Conclusions 
follows in section (l6l). 



2. Data set and numerical simulations 

As briefly highlight in the introduction, our main purpose is to investigate intermit- 
tency of the Lagrangian/Eulerian velocity field in the inertial range for a weakly com- 
pressible isotropic turbulent flow (see figure ([T]) for a colorful rendering). The numerical 
simulation have been performed using the FLASH code developed by the Flash Center 
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at the University of Chicago ( Fisher et alJ ( 2007 )). The equation of motion are: 

'dtp+V ■iup) = 
dtipu) + V • (uup) = -VP + F 
dt{pE) + V -[uipE + P)] =0 

P=i^-l)pU 
E = p{U + ^u'^) 



(2.1) 



where p is the mass density, u the velocity, P the pressure, E the total energy density, 
U the specific internal energy, and 7 is the ratio of the specific heats in the system. 
Finally, the fourth equation in (|2.ip is the equation of state for our system. The effect 
of the large scale forcing F gives rise to a turbulent flow whose energy is transferred 
from scale Lq towards small scales. The energy input ^ (Pxvl ■ F produces an increase 
of the internal energy U , which grows in time. One can easily show that the quantity 
J dx^PdiVi represents the energy transfer from kinetic to internal energy of the flow. 
The sound speed increases in time as well, while the average Mach number is order 0.3. 
The numerical simulation is done for isotropic and homogeneous forcing with a resolution 
1856"^. The integration in time have been done for three eddy turnove r times after atran- 
sient evolution. More details on the numerical code can be found in lBenzi et al.l (j2008l ). 



The code develops density gradients due to compressibility. Despite compres sible effects 



which l ead to highly non-t rivial scaling of density and entropy fluctuations (jBenzi et al 



( 20081 ) : IPorter et all (l2002l )). i;hc net feedback on the velocity scaling is very weak, if any. 
Here we show that to the best of our testing ability, velocity fluctuations in the inertial 
range are indistinguishable from those measured on incompressible fluid turbulence. 
Although the integration is formally inviscid, there is a net energy transfer from the 
turbulent kinetic energy \/2pu^ to the internal energy. Thus we may think that an ef- 
fective viscosity is acting on the s ystern . In order to estimate it, we proceed as if the 
Kolmogorov 4/5 equation (see lFrisch (1993)), valid for incompressible turbulence applies 
to our case -with effective parameters: 



^cff^ 



, d (2) 

'eff^^L I 



A fit of our data with this formula gives, 

,3 / 1/4n 



r). (2.2) 



8.3-10 ^ , which corre- 



0.054, 

sponds to a Kolmogorov scale, 77 = (^cff/^eff )' equivalent to roughly half grid cell and 
to R\ ^ 600. The dynamical effects of the effective viscosity is however different from 
what one usually observes in the Navier Stokes equations, i.e. the dissipation range does 
not behave similarly to the Navier Stokes solutions. Thus, while on the energy flux (i.e. 
the third order structure functions) we still observe an effective dissipation, the dissipa- 
tion range is changed according to the specific mechanism employed in the simulation, 
namely the increase of inter nal energy and the artificial numerical viscosity added up to 



smooth too high gradients (jFisher et al.l (j2007t )). Let us quantify the importance of the 
small compressibility on the statistics of the velocity field. For isotropic and incompress- 
ible fields, there exists an exa ct relat i on wh ich connects second order longitudinal and 



transverse structure functions ( Frisch (|l995l) ) 



r d (2) 



(2.3) 



this equation is useful because any deviations from it gives a quantitative hint on the 
cumulative importance of anisotropy and compressibility, scale-by-scale. We show in the 
left panel of Fig. ^ the comparison between Sj. ' (r) and its reconstruction via the RHS 
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Figure 1. Example of Eulerian and Lagrangian rendering. Both the intensity of the Eulerian 
velocity field at a given time and the Lagrangian evolution of a bunch of particles with tra- 
jectories ending at the time of the Eulerian snapshot are shown. Notice that the Lagrangian 
particles have an initial smooth evolution because we also show the initial transient time when 
the underlying Eulerian field was chosen smooth and with low energy. The bunch of particles 
was chosen such has to encounter a vortex filaments during their evolution. Figure courtesy of 
B. Gallagher. 




Figure 2. Left: Test of the isotropic and incompressible constraint, we show the LHS and 
the RHS of (|2.3p . In the inset we have a percentage estimate of the breaking of the relation: 
(RHS-LHS) /LHS. Notice that percentage- wise the relation is well verified, within 5% of accuracy. 

Right: third order longitudinal structure functions, S'jf'(r) with the power law exact relation 

coming from the 4/5 law (|Frischl (|l995f )). Inset: compensated plot, Sf\r)/r. 



of (|2.3p . As one can see the agreement is good. Combined effects due to anisotropy and 
compressibility are less than a few percent for second-order statistics (see also Fig|3]for 
higher-order statistics). 
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3. Intermittency and anomalous scaling in the Eulerian velocity field 

We start our analysis by measuring the scaling behavior of the Eulerian structure func- 
tions. In the right panel of Fig. ^ we show the third order longitudinal ESF as measured 
in our DNS. As one can see, in log-log plot there is a clean scaling range over more than 
one decade extending from r/ij € [10 : 100], with a slope close to the scaling prediction, 
^^'(r) ~ r. In Fig. © we show a typical log-log plot of Longitudinal ESF for p — 2,6, 10. 
Again, a qualitative scaling behaviour is clearly detectable. The high resolution and the 
highly isotropic forcing used in our DNS, allows also for a precise quantitative assessment 
of the importance of statistical and anisotropy effects on the averaged quantified. In the 
right panel of the same figure, we show the estimate of the statistical fluctuations (top 
panel) by plotting the ration between two Structure functions averaged either over the 
whole statistics or over one half of it. As one can see, statistical fluctuations become at 
maximum of order 10 — 20% at small scales (where intermittency is more severe) and 
only for high order moments p = 6 and larger. Inside the inertial range, r > IO77, they 
are order 5% or even less. Similarly, in the bottom panel we give a quantitative estimate 
of anisotropic residual effects, by plotting the ratio between longitudinal measurements 
over two different orthogonal directions. Here, as expected, deviations are larger close 
to the integral scale, reaching a maximum of 10% for high order moments. As a result, 
the combined effect of statistical and anisotropic fluctuations for the Eulerian scaling 
are small and they result in an estimate of the error bars that is within the size of the 
symbols for the log- log plot showed on the left panel of Fig. ^ . 

One of the pluses of high resolution DNS is the possibility to go beyond log-log plot, 
analysing scaling properties locally, with high precision. In order to do that, we analyse 
scaling by using the magnifying glass of local scaling exponent, (LSE), given by the 
log-derivative of Structure Functions: 

(,) dlogSg^. (P) dlog4^^(r) 
^"^^ - dlogr ' - dlogr ■ 

The advantage to rely on LSE is twofold. First, it allows for assessing statistical prop- 
erties locally, removing large-order non-universal contributions coming from the overall 
prefactors in the SF. Second, and more importantly, they assess scaling without the need 
for any fitting; in other words, they are just the outcome of a measure and they do not 
depend on an arbitrary definition as the extension of the inertial range, in order to de- 
fine the degree of intermittency. For instance, let us consider, the scaling of fourth order 
longitudinal Flatness. It is easy to see that it can be rewritten, independently on any 
scaling assumption and for all separations r, as: 

In other words, the presence of intermittency, i.e. a non trivial scale-dependent Flatness 

behaviour, is directly measure by how much the ratio of LSE, Cl'''(^)/CF^(^) different 
from its dimensional value, 2, scale-by-scale. Let us start to analyze the LSE for the 
Longitudinal and Transverse ESF. In figure Q we show the LSE for p = 2, 4, 6, 8 and in 
figure (O we show the LSE for p = 10. Let us first notice that Longitudinal and Transverse 
statistics have significantly different viscous cut-off, with transverse fluctuations having 
the tendency to have le ss important viscous damping with respect to the longitudinal 
one ( Gotoh et al This is an effect already present for low order p = 2 (bottom 



left panel) and therefore certainly induced, at least partially, by geometrical constraints 
as the ones discussed about eqn. (|2.3|1 . Beside this remark, let us also notice the small 
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Figure 3. Left: log-log plot of Longitudinal Eulerian Structure Functions, S^j^\r) forp = 2, 6, 10. 
Right: estimate of the anisotropic fluctuations, ([^^(x + r^) — Ui(x)]*') / {[uy{x + ry) — Uy{-x.yf), 
for p = 2, 6, 10 (bottom). Estimate of the statistical fluctuations, obtained by plotting the ratio 
of S^l\r) calculated over the whole statistics or over one half of it (top). The sum of the two 
error sources leads to error bars that are of the order of the symbol size on the left panel. 



-in amplitude- but important -in principle- difference in the scaling exponents form 
longitudinal and transverse. Up to p = 4 the LSE are almost coinciding, within error bars, 
in the range of scales r/ry G [20 : 80]. For larger orders, we start to detect a difference, 
with the transverse being systematically below the longitudinal. This difference becomes 
significantly visible for p = 10 as shown separately in Fig. ([5]). Note that this statements 
is not the outcome of any fitting procedure and what is shown in figure (|4]) and (O 
is a measurements independent on any theoretical interpretation. To decide about the 
existence -and extension- of a power-law scaling range needs some fitting. The value of 
the global scaling exponent, Cl'^ Ct''' would then be given by the average of the LSE 
in the region where they are close to constant. An error bar on the mean value can 
be then estimated on the basis of the oscillation induced by statistical and anisotropic 
effects in the fitting range and by the change as a function of the extension of the scaling 
range used to make the fit. Such a fitting procedure (details in the caption) leads to 
the summary for the global Eulerian scaling exponents depicted in the right panel of 
Fig. ([5]) and summarized also in table 1, where we compare our results also with other 
DNS obtained with fully incompressible Navier-Stokes equations at comparable Reynolds 
numbers. From a theoretical point of view, we know that there are exact constraints that 
fix ^ = Ct ^ for p = 2 (for p = 3 a similar constraint fix the scaling of third order 
longitudinal with mixed second order transverse and first order longitudinal averages). 
Nevertheless, there is not any strong theoretical argument suggesting the possibility that 
for p > A one should expect different scaling. So, the result shown in Figs. (jUlS]) is 
not fully understood. Whether the difference between longitudinal and transverse ESF 
will shrink going to higher Reynolds number remains to be investigated, and it is an 
important open question. The good news is that the value found in our numeric are 
in perfect agreement with previous results (see table 1), making us confident that they 
are robust and not dependent on compressibility. Even the relative scaling behaviour, 
obtained by plotting each structure functions versus the thi rd order one , a procedur e 
known as Extended Self Similarity (ESS) in the literature ( Benzi et al. ( 1993 . 1996f )) 



does not change much the above picture. In table 1 we give also the results of the LSE 
for both Longitudinal and Transverse ESF estimated by using ESS. Again a small spread 
between longitudinal and transverse is measured. Nevertheless, it is interesting to notice 
that while the ESS procedure applied to longitudinal statistics reduces considerably the 
error bars with respect the usual LSE, for the transverse one there is not a clear gain 
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1 10 r/ri 100 1 10 f'"! 100 




1 10 r/ri 100 1 10 r/rj 100 



Figure 4. Local scaling exponents, Cl (r), Ct of Eulerian Lagrangian and Transverse struc- 
ture functions for p = 2,4, 6, 8. Error bars are obtained by estimating the residual anisotropy, 
and they are given by the spread between measurements on different directions, r/r — x,y,z. 
Statistical error bars are smaller. The horizontal straight line corresponds to the k4f dimensional 
prediction, (^'''^ = p/3. 

in adopting ESS. This may suggest the possibihty that the different scaling behaviour 
detected between, S^j^\r) and s!^\r) maybe due to different sub-leading corrections 
to the main leading scahng, contributing more to the transverse scaling than to the 
longitudinal. Due to the lack of any hints on the possible sub-leading correction we refrain 
here from entering in to a fitting procedure with too many free parameters. The issue 
whether the different in the LSE shown in Fig. ^ is the signature of a true mismatch 
between the scaling properties of longitudinal and transverse fluctuations or the result of 
a superposition of leading and sub-leading power laws with coinciding exponents between 
longitudinal and transverse but with different prefactors, remains open. 



4. Intemittency and anomalous scaling in the Lagrangian velocity 
field 

We now turn our attention to Lagrangian particles. As it is known, LSF shows very 
strong intermittency. In figure ([6]), we show the probability density function of drVi, for 
different values of r and averaging over the different velocity components: for small r the 
pr obability density exhib its stretched exponential tails, as also measured experimentally 
in ( Mordant et al. I (l200lh l The huge statistics at our hands, up to 16 millions trajectories 
of length up to one large scale eddy turn over time, Tj,, allows for detecting fluctuations 
up to 40 standard deviations, for the highest frequency. Concerning the LSF, we show in 
the same figure (right panel) a first overview on a log-log scale for moments up to p = 8. 
From this panel we can already extract some conclusions. First, scaling in Lagrangian 
framework is not as good as the Eulerian one, as one can judge by the naked eye. This 
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Left: longitudinal and transverse local scaling exponents, the same of fig. Q 
= 10. Right; summary of the Eulerian scaling from table 1. The two solid lines 
represent a fit to the scaling exponents with two different Mult i Fract a l spectrum DlI Ii) 

' ■ T)) 



and Drih) both having a log-Poisson statistics, (jShe &: Levegud (|1994l ): iDubrull^ (|199 

DL,T{h) 



3(h-h') 



log 



3{h*-h) 



1 + 3 — T . The free parameters are chosen such as 



to minimize the fitting differences between Longitudinal and Transverse. The only parameter we 
changed is the intermittent factor = 0.6 and /3t ~ 0.4 while we kept the highest singularity 
to be the same for the two cases, h* = 1/9. The fractal dimension of the h* is then fixed by 

the requirement that ("j^''' = ("y^' = 1 and must be chosen dl ^ = (1 — 3h*)/{l — Pl,t)- Another 
alternative approach for the fit would be to fix (3 and change the value of the highest singularity 
between longitudinal and transverse fluctuations. 



Ap) 

ST 



Ap) 



[G02] 



Ap) 
ST 



[G02] 



2 0.71 ± 0.02 0.71 ± 0.02 0.70 ± 0.01 0.71 ± 0.01 0.69 ± 0.005 0.71 ± 0.01 

4 1.29 ± 0.03 1.27 ± 0.05 1.29 ± 0.03 1.26 ± 0.02 1.28 ± 0.01 1.26 ± 0.01 

6 1.78 ± 0.04 1.68 ± 0.06 1.77 ± 0.04 1.67 ± 0.04 1.75 ± 0.01 1.68 ± 0.03 

8 2.18 ± 0.05 1.92 ± 0.10 2.17 ± 0.07 1.93 ± 0.09 2.17 ±0.03 1.98 ± 0.1 

10 2.50 ± 0.06 2.10 ± 0.20 2.53 ± 0.09 2.08 ± 0.18 2.5 ± 0.05 2.25 ± 0.15 

Table 1. Estimate of the global scaling exponents, out of the figures for the local scaling 
exponents The first and s econd columns refer to our DNS, the third and fourth to the 

data published in (jGotoh et al.l (|2002l )) from a DNS at comparable Reynolds number. The last 
two columns correspond to the ESS estimate of our data, using the third order longitudinal ESF 

ip) 

as a reference for the scaling of Si and using the third order transverse ESF for the scaling 

of S'^^(r). Error bars are obtained by summing the uncertainty obtained from estimating the 
scaling in three different spatial direction (anisotropy contributions) and by changing the scaling 
range in the interval r £ [10 : 100]?;, where the global exponent is evaluated. 



is a common feature of all Lagrangian statistics, and was already observed in many 
other experimental and numerical previous works. The simplest explanation is that the 
finite Reynolds effects are more important in Lagrangian domain: dimensional estimate 

1/2 

gives for the Lagrangian inertial range extension the sealing Tl/t,, ~ , while in 

3/4 

the Eulerian case we have L/rj ^ R-^ . In the inset of the same figure, we show the 
importance of anisotropic fiuctuations, by comparing the ration between two LSF on two 
different components, s'i'^ (r) /S^^ (r) for p = 4, 8. As one can see, the importance of 
anisotropic fiuctuations is again of the order of 10% at most, for the highest moments. 
Such correction are of the order of the symbol size of the LSF averaged over the three 
components shown in the left panel. The poorer scaling in the Lagrangian domain, does 
not allow for a systematic assessment on the local scaling exponents as done for the 
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Figure 6. Left: Lagrangian PDF for a single component velocity increments along particle 
trajectories, over different time increments, r £ [2 : 400]r,,. Right: log-log plot of LSF for 
p = 2,4, 6, 8. Inset: estimate of the anisotropic statistical degree: ratio between two different 
components of LSF, {{S-rVxY) / {{SrVyY) , for p = 4 (circles) and p — S (triangles) 



Eulerian field. Here, if one wishes to retain good local properties, it is necessary to resort 
to the ESS method, plotting the local scaling exponents relative to one moment versus 
a reference one (here taken the second order): 

which are related, in presence of pure power law scaling, to the scaling exponents of the 
LSF defined from (|1.2p by the obvious relation: 

X^^\t) const. =C^P^/C^^l (4.2) 

Let us stress nevertheless, that the importance of the local exponents (|4.ip goes much 
beyond their interpretation as a simple proxy of the ratio between the structure functions 
exponents (|4.2p . Indeed, they give a clear and simple way to assess the importance of 
intermittency in the Lagrangian domain being the only quantities entering in the scaling 
properties of Lagrangian hyper-Flatness: 

^'"Hr) = = (5(^)(,))x-)(.)-./2. (4.3) 

Again here the same comment made for the Eulerian case is in order: via the Lagrangian 
hyper-Flatness (|4.3p we are able to assess intermittency in a quantitative way, free of any 
fitting ambiguity, without having to assume power law properties, by simply checking 
the difference between x''^^''(''') ^-^^ p/2, scale-by-scale. In Fig. ([7]) we show the hyper- 
Flatness local exponents for p ~ 4,6,8,10. The scaling range now is good, for r/r,, g 
[20 : 200]. As one can see, there are two remarkable facts. First, Lagrangian scaling is 
much more intermittent then the Eulerian one, if measured as a deviations from the 
dimensional scaling, xip) =p/2. For example, for p = 10, we measure x'^^^'^ = 2.45 ±0.35 
which is 50% off from the dimensional value x^'^^^ = 5. The second fact is the strong 
dip measured across the v iscous scale, for time lags t/t,, G [1 : 10]. This is similar to 
the Eulerial "bot tleneck" dLohsd (Il994l)), w hile for Lagrangian structure functions was 
observed first in ( Mazzitelli fc Lohs d ("2004^) and interp r eted as induced b y the presence 



observed hrst m [|Mazzitem &: Lohs e [ z!UU4')) and mterp r eted as mduced b y the presence 
of small scale vortex filaments in [Biferale et~all |2005l) : iBec et al.l (|2006l ). The value of 



the Lagrangian exponents measured in our DNS ar reported in table (2) and compared 
with previous measurements on other DNS or experiments. Let us stress that there are 
not other measurements with comparable statistics and Reynolds number, as the one here 
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Figure 7. Lagrangian Local scaling exponents for hyper-Flatness, x'"'''' (''')) for P = 4,6,8,10. 
Notice the extremely good scaling behaviour for low moments, which deteriorates only mildly 
for p = 6, 8. Notice also the dip region close to the viscous t i me scales, an eff ect interpreted in 
terms of trapping into vortex filaments ([Jimenez et alj (| 19931 ) : lBec et all (|2006f )). The horizontal 
straight lines correspond to the dimensional scaling x'''^ = p/2. 



X 



(p) 



10 



our DNS {Rx ~ 600) 1.66 ± 0.02 2.10 ± 0.10 

EXPl {Rx ~ 800) 1.47 ± 0.18 1.73 ± 0.25 

EXP2 {Rx e [500 : 1000]) 1.56 ± 0.06 1.80 ± 0.20 

DNSl {Rx ~ 320) 1.51 ± 0.04 1.76 ± 0.11 

DNS2 {Rx ~ 75) 1.56 ± 0.03 1.82 ± 0.08 



2.33 ± 0.17 
1.92 ± 0.32 



2.45 ± 0.35 
1.98 ± 0.38 



1.92 ± 0.14 1.93 ± 0.3 



Table 2. First raw: summary of the local scaling exponents for Lagrangian statistics as extracted 
from our data shown in Fig. (IT}. Error bars are evaluated out of the anisotropic statistical degree, 
estimating the LSE out of the three different velocity components and by changing the fit range 
in the interval r G [20 : 200]. The other rows correspond to other data available for high 
order moments from both experiment s at comparable Reynolds numbers and numeric al works 
at smaller Reynolds numbers. EXPl (IXu e t al.' f200^); EXP2 ([Mordant et all (|2006l )): DNSl 
([Homann ct al. (2003)); DNS2 ([Mordant et al. (2006)). Notice how some of the previous data 
slightly underestimate the value of the scaling exponents for p = 6, 8, 10. This is due to the 
fact that, due to the limitations either in the statistics or in the Reynolds number, it was not 
possible to assess scaling properti es scale- by-scale as don e here . For example, the fit leading 
to the numbers in IXu et aP ([20061 1 and in I Homann et al.l ([20071 1 (second and fourth raw) was 
done in the range r £ [3 : 6]t„, well insi de the dip region. Similarly, the Reyolds number of the 
numerical data in ([Mordant et al.l ([2006l 'l) (fifth row) is very low, probably the strong saturation 
in the values for p = 6, 8, 10 reported there is due again to the presence of the dip region spoiling 
the scaling in the short inertial range present at Rx ^ 75. 
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presented for high order moments. A global fit in log-log of the scaling properties does 
not allow to disentangle the inertial range from viscous-inertial contribution (affected by 
the dip) as they emerge from fig. ([7]). This explain the slightly systematic underestimate 
of the scaling exponents from all the previous studies with respect to our measurements. 
The only da ta presented pre viously for p = 8, 10 at comparable Reynolds are those 
published in (|Xu et all (l2006l) ) where the scaling properties where estim ated very close 
if not inside, the dip region, r £ [3 : 6]tj,. Similarly was d one for the data in Hom ann et al 



([20071). The numerical data for p = 8, 10 presented in (jMordant et al.. (.2006)) are also 
understimating the value of the scaling properties with respect to our data. Probably, 
in this case, the very low value of the Reynolds number {Rx 75) is such that the 
whole scaling region is dominated by the viscous dip. Considering all these pitfalls, we 
may conclude that when scaling properties are analysed in the correct range and where 
statistical fluctuations are not too high, there exis ts a good uniyersalit y in Lag rangian 
statist ics, as also tested quantitatively for p = 4, 6 in Biferale et alj ( 20081 ) : lArneodo et all 
(|2008^ . As for universality of high order moments, from p = 8 and up, we need to wait 
for other data with statistical properties and Reynolds number comparable with those 
of this study, to decide. 



5. Multifractal: a link between Eulerian and Lagrangian statistics 

Multifractals have been introduce 25 years ago to explain deviat i ons to the K41 



scaling for Eu l erian isotropic turbulent fluctuations (jParisi fc FrischI (jl985f ). see also 



Boffetta et al.l (j2008l) for a recent r eview). Since t h ose works, they ha ve been extended 



to describe also vel ocity gradients (Nelkin ( 1990l): Benzi et al. ( 1991 )): multi-scale ve- 



locity correlations ( Belinic her et al.' ('l998'l: 'Benzi et al.l (Il998l)') and fluctuations of the 



Kolmogorov viscous scale ([P aladin & Vulpiani (1987); Menevcau (.1996.' ): .Biferal e ( 20081 ): 
iFrisch fc Vergassola (|l99ll) ). Let us notice that other attem pts have been mad e to in - 
troduce fluctuations o f the Kolmogorov scale in turbulence (jLVov fc Procaccial ( 19961 ): 
Yakhot fc Sreenivasa 3 (|2004. 2005)), most of the m are equivalent or gi ve anyhow results 



almos t undistinguishable from the multifractal approach followed here (jBenzi fc Biferale 
(|2009[ )). 

In principle, one is free to develop different, uncorrelated, m ultifractal description for 



Eulerian and Lagrangian quantities ( Chevillard et al. ( 20031 )). On the other hand 



reality, Eulerian and Lagrangian measurements are of course intimately linked. One is 
therefore naturally tempted to develop a unified Multifra c tal description, able to describe 
both temporal and spatial fluctuations ( Borgad ( 2003f l: Boffetta et al.l ( 2002 )). In the 
following we describe such development starting from the inertial range and finally for 
the viscous range. 

5.1. Inertial range 

Concerning Eulerian MF in the inertial range, the idea goes back to lParisi fc FrischI (|198,^ 
and reads as follows. Suppose we have different velocity fluctuations with different local 



Holder exponents, Sr 



for scale separation in the inertial range. Suppose that the 



set where the velocity field has a /i-exponent is a fractal set with dimension D{h). Then 
the ESF can be easily rewritten as an ensemble average over all possible h fluctuations: 



{{5rUY) 



dhr^^r^-^^^^ 



(5.1) 



where the factor r^ gives the probability to fall on the fractal set with ^.-exponent at 
scale r and where we have neglected, for the moment, any difference between longitudinal 
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and transverse, calling the generic Eulerian increment, 6rU. The last passage in (|5.ip is 
obtained in the saddle point limit r — > 0: 

C'^P^ = min {hp + 3 - D{h)) . (5.2) 

h 

It is easy to imagine many different possible D{h) distribution -built in terms of random 
cascade- leading to a set of exponents C'^-' close to those measured. In previous section 
we have shown that both our DNS and previous DNS suggest the possibility to have 
different scaling exponents for Longitudinal or Transverse SF. Therefore we must allow 
to have two different set of fractal dimensions, DL{h) and Dxih) describing the statistics. 
For example, in the right panel of Fig. ([5]) we show the result of the saddle point (|5.2p 
obtained with two different fractal spectra, tuned to fit the empirical scaling exponents 
as found from our DNS (see caption for details). 

One po ssible simple w ay to link the Eulerian and Lagrangian statistics is to follows the 



ideas in iBorgad (|2003f ). which we summarize hereafter. As far as scaling is concerned, in 
3d turbulence we may imagine that the only relevant time in the inertial range is the 
local eddy turn over time: 

r(r) ~r/(5,u) -r^-''. (5.3) 

Then, we may assume that the Lagrangian velocity increment over a time lag r must be 
of the order of the corresponding Eulerian velocity increments over a scale r with r and 
T connected by (|5.3p : 

SrV '^5rU T = r^~'\ (5.4) 

By using only this simple statement one is able to link now the LSF to the ESF, by 
simply writing: 

f hp 3-0(h) ^(p) 

{{Srv)P) = / dhT—r^^ , (5.5) 

Jh 

with now the Lagrangian exponents given by the saddle point estimate: 

t(p) fhp + 3-D{h) 



It is important to notice that via the translation factor (|5.3p we obtain two different sets 
of Lagrangian and Eulerian exponents but given by the same D(h) curve. The exponents 
are different because the local eddy-turn-over is itself a fluctuation quantity, depending 
on the local scaling exponent, h. 

Once given the _D(/i), extracted from a fitting of the Eulerian statistics, it is tempt- 
ing to use it as an input in (|5.6p to get a prediction with no free parameters for the 
Lagrangian scaling. Indeed, because the Eulerian statistics is different depending if one 
takes longitudinal or transverse fluctuations, and because it is very natural to think that 
during the Lagrangian evolution both fluctuations are felt, one may imagine to get two 
different prediction for the Lagrangian exponents: 

^(p) - ..i. f hp + 3-DL{h) \ . , _^,^fhp + '3- Drih) 



""t \ i—h ) ' ^ = Y—h ) ■ (^-^^ 

depending which Eulerian statistics one uses, DL{h) or Drih). In fig. ([5]) we show the 
comparison between the empirical Lagrangian exponents as extracted from our DNS (see 
table 2) and Fig. ([7]) and the two prediction one get using the first or second expression 
in (|5.7p . As one can see the agreement is good, indicating that the simple, but non 
trivial, bridge relation (|5.3p is working well and capturing the main scaling properties 
of both Eulerian and Lagrangian statistics. From figure ([5|) one may see that only for 
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Figure 8. Summary of the Lagrangian Scaling exponents as measured in our DNS in the inertial 
range of time scales (circles), together with the prediction obtained from the Eulerian statistics 
by using the bridge relation either with the longitudinal Eulerian statistics (upper limit of the 
shadowed area) or with the transverse one (lower limit of the shadowed area). 



the highest order p = 10 one start to see a trend deviating from what predicted by the 
Eulerian-Lagrangian bridge relation, although within error bars it still holds. If this is 
the signature that very intense Lagrangian fluctuations bring new information beyond 
the ones collected by (j5.3p is an open and interesting question, which is not possible to 
address till new empirical and numerical data will allow to access higher order statistics 
for both Eulerian and Lagrangian domains. 

5.2. Dissipative Effects 

The possibility to assess Lagrangian data, both numerically and experimentally, opened 
the way to investigate also dissipative and sub-dissipative time scales. For Eulerian 
statistics is very difficult to have reliable data at scales smaller then ry, due to ei- 
ther limitati ons induced by th e probe size, for experime n ts, or by the numeri c al res - 
olution (see ISchumached (|2007^ : ISchumacher et all (|2007l ): IWatanabe fc Gotohl (|2007t ): 
Yamazaki et al.l (j2002[ ) for recent DNS meant to address this issue). For Lagrangian 
quantities , the development of non-intr usive experimental techniques (mainly based on 
fast CCD iToschi fc Bodenschat d (l2009l) ) and the advancements of numerical algorithm 
to track particles, allowed for assessing time fluctuations of the order of O.lr^ even at 
high Reynolds numbers. Now, new question about the statistical fluctuations for r ~ 
or smaller, arise. Recently, a lot of attention has been payed on the dip region visi- 
ble in the local Lagrangian s caling exponents as s hown in fig. ([7]) both from numeri- 



cal and experimental studies ( Arneodo et al. ( 20081 )). In the last twenty years, a lot of 
work has been done to i nclude viscous fluctuations i n the multifrae t al theory - i nitiall y 



for Eulerian statistics- IPaladin fc Vulpianil (fl987h : iNelkin I (|l990f ): iMeneveaul (|l996h : 
Frisch fc Vergassolal ( 199ll l: IBenzi et al.l ( 199ll ). and more recently also for Lagrangian 
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Biferale et alJ ( 20041 )). A simple way to accomplish 



this is to take into account of the fact that, according to the MF theory, the viscous 
scale is not a fixed homogeneous quantity, but fluctuates wildly due to intermittency. 
Moreover, viscous and inertial fluctuations are linked by the requirement that viscous 
effects happens wh e n the local Reynolds number becomes 0(1) as first remarked by 
Paladin fc Vulpianil (|l987[) : 

5riU rj 



0{1) 



V 



1+h 



(5.8) 



where we have taken for simplicity, L = l,Uo — 1 and v = Re~^ . Using the bridge 
relation (|5.3p we may immediatel y find out also the e quivalent fluctuations relation for 
the Kolmogorov dissipative time (jBiferale et al. 



l-h 



(I2004D ): 



(5.9) 



Taking in t o account this fact, one may use a Batchelor-Menev e au like parametrisatio n 
( Meneveaul |l996l ); ISirovich et all (|l994h ; IChevillard et all (jioosl ) : I Arneodo et all (|2008l )) 
of velocity increments for all time lags, including the fluctuating viscous time in the game: 



^rV = Vo — 



Tr 



(5.10) 



where a being a free parameter controlling the crossover around t^t^j, and Vq the root 
mean square large scale velocity. It is easy to realize that the above expression has the 
correct asymptotic behaviour for very small time lags StV ~ to, where a = Sr^vTrj is 
the acceleration; and in the inertial range. S^-v ~ r''/(^~''\ as described in details in 
Arneodo et al. Iliooi). In order to get a prediction for the behaviour of the LSF we need 
also to introduce dissipative effects in the MF probability. This can be done as follows: 



Tl 



Tl 



(5.11) 



where 2 is a normalizing function and D{h) the fractal dimension of the support of 
the exponents h (either DL{h) or Dxih) if one wants to link the Lagrangian dissipative 
physics to the Eulerian one). At this point, given the Reynolds number, we are left with 
two parameters {a and a multiplicative constant in the definition of r,,). In order to 
compute the LSF at all time lags one needs to integrate the expression (|5.10|) with the 
weight given by (|5.1ip over all possible /i-values: 



iiSrvr) 



dhPh{T)[SrVih)]P. 



(5.12) 



Clearly, the above picture goes back to the inertial values when we have t/Tl —>■ and 
T ^ Tri, i.e. when we can neglect viscous effects and we can perform the saddle point 
estimate. Here we do not want to do that, keeping all corrections, and compare the simple 
parametrisation (|5.10p . dressed with the MF theory, with our data. First, let us see what 
are the qualitative results obtained from (|5.10p at changing the free parameters on it. 
In fig. ^ we show the results for x^p\t) as obtained from the global MF description 
(|5.12p for p = 4, 6, 8 with a fixed D{h) at changing the free parameter a = 4 (left panel) 
or a = 2 (right panel). From this figure we learn two things. The MF description is very 
close, already qualitatively speaking, to the real data -compare with panel corresponding 
to p = 4 of Fig. d?]). In particular, the dip present in the data is also captured by the 
MF. Moreover, at changing the free parameter a one may deplete or enhance the dip 
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intensity, as can be see comparing left and right panel of Fig. ([9]). Let us also notice that 
the dip region is entirely due to the presence of the fluctuating viscous cut-off, t,j(^), as 
can be seen in the inset of right panel of Fig. Q , where the same data but with a fixed 
Trj are plotted and the dip has disappeared. 

Second, from fig. ([9|) one sees that for any finite Reynolds numbers, there are corrections 
to scaling induced by the large time lags T^; as a consequence, in the inertial range the 
local exponents are not exactly constant. In order to give a qualitative assessment of the 
importance of this finite Reynolds corrections as described by formula (|5.10p we plot in 
the left panel of fig. (fTU|) the values of x^'*-' (''■) for three different Reynolds numbers. As one 
can see, at increasing Reynolds, the inertial range values becomes flatter and flatter, as 
they must do. At difference from what happens in the infinite Reynolds number limits, for 
any finite Reynolds, the whole support of D{h) will play a role in the integration ()5.12|) : 
for not too small t/Tl the saddle point estimate (|5.6p becomes less and less accurate. 
As a result, the whole shape of the D(h) curve becomes important for any order. 
Indeed, the typical D{h) curve from any MF theory is a convex function like the one used 
to fit the Eulerian Longitudinal exponents depicted in the inset of right panel of fig. (fTO]) . 
The left part of D{h) curve, with h £ \hmin ^ hnfar t.-], is connected to the positive order 
of the Eulerian Structure Functions I FrischI ( 19951 )'). The right side, h G [hpeack, hmax], 
is connected to negative moments. Negative moments are generally hill defined and on e 
needs to use inverse statistics to assess the right part of the D{h) curve (j Jensen (Il999l) : 
Biferale et al. 1 (119991) 1. Because of that, we have a good control only the left side, where 



most of the previous experimental and numerical studies focused. The extension of the 
integration in ()5.12p over the h-value on the right of the peack of D{h) is therefore not 
based on solid experimental or numerical basis. The extension to the right part of a given 
functional form for D{h) (log-Poisson, log-Normal etc. ), obtained by fitting only positive 
moments of Structure functions, is arbitrary. 

In the right panel we show the dependence of the final expression (|5.12p from the 
extension of the integration over the h range. As one can sec, the inclusion of ft,-value 
falling on the right of the peack gives an important contributions to the behaviour of 
the local slopes for t ^ Tl, confirming the importance of the whole D{h) if one want 
to control the finite-Reynolds corrections to the saddle point estimate (|5.6p . For a more 
detailed investigation of the property of the MF formulation (|5.12[) also compared with 
other DNS a nd experimental results, the reader can consult the recent collection of data 
pubhshed in (jArneodo et all (|2008h ). 



6. Conclusions 

We have reported about Eulerian and Lagrangian statistics from high resolution Direct 
Numerical Simulations of isotropic weakly compressible turbulence. Reynolds number at 
the Taylor microscale is estimated to be around 600. The high number of particle tra- 
jectories explored, together with the high resolution Eulerian flow allowed us for the 
first time to study Lagrangian intermittency up to moments of order p ~ 10, at large 
Reynolds numbers. The four main results we have presented are: (i) Eulerian Longitudi- 
nal and Transverse statistics have slightly different scaling properties, in agreement with 
previous incompressible DNS data at comparable resolution, (ii) Lagrangian statistics 
is more intermittent than Eulerian ones, percentage-wise with respect to the observed 
deviations from dimensional scaling and its quantitative value can be captured by using 
a simple multifractal bridge relations, (iii) In order to correctly asses Lagrangian scaling 
properties is mandatory to make a scale-by-scale analysis, in order to disentangle the 
dip region close to viscous scales from the inerital range interval. A lack of such anal- 
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Figure 9. Results for the local Lagrangian exponents, x'^'H''') ^-^ obtained from the finite 
Reynolds multifractal formula (|5.10I5.1H5.12|I for a given Reynolds number and p = 4,6, 8, 
at changing the free parameter, a, that sets the importance of the dip region. Left; a = 4; 
Right: Q = 2. In the inset of the right panel we also show, x'^' (''") obtained from the same 
multifractal formula but keeping fixed the dissipative time scales, r^(ft) = const., in this latter 
case the dip region disappears. 
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Figure 10. Left panel: Reynolds effects on the local exponents predicted from the MF in- 
terpolation formula (|5.10I5.1H5.12|) . we kept fixed everything except the Reynolds numbers. 
Right: sensitivity to the whole shape of the D{h) curve of the x'''' (''") as predicted from the 
MF interpolation formula. We show two cases. One case is obtained integrating only on fluc- 
tuations in the left region of the D{h) curve with respect to its peack, h G [hmin, hmax], with 
hmax = hpeack = 0.38, the secoud case consider also smoother fluctuations up to hmax ~ 0.5. The 
typical log-Poisson shape of the D{h) used is depicted in the inset. 



ysis may lead to an underestimate of the Lagrangian scaling exponents, because of the 
strong intermittency of the dip interval, dominated by viscous effects, (iv) Multifractal 
phenomenology, is able to describe in a unified way: dissipa tive effects, incrtial range 
fluctuations and large scale corrections (jArneodo et al ] (I2OO8I) ). 

This study leaves open a few important questions. First, we need to wait for higher 
Reynolds numbers data, in highly isotropic ensemble, to resolve the riddle about Lon- 
gitudinal versus Transverse Eulerian fluctuations. Second, we need to collect also data 
with higher Lagrangian statistics to understand whether the Eulerian-Lagrangian bridge 
relation remains valid for higher moments. Third, it is important to improve experimen- 
tal and numerical accuracy to measure small-scales and small-time fluctuations, where 
highly non trivial physics is developing as shown by the strong enhancement of local inter- 
mittency in the dip region r G [1 : lOjr,,. Such strong enhancement of fluctuations around 
the viscous scale is due to local fluctuations of the dissipative cut-off, which reflects in to 
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the existe nce of different viscous effe c ts for different moments and different correlation 
functi o ns (IFrisch fc Vergassola (jl99ll ): iL'vov fc Procacci 1 (HoOfih i lYaThot fc Sreenivaii[^ 



(|2004[ ): iBenzi fc Biferald (|2009r )). ft is well described, within the multifractal theory, by 
the Paladin- Vulpiani relation (|5.8p . as shown by the fact that the dip region disappear by 
removing the viscous fluctuations in (|5. 10115. lT|) . see inset of the right side of fig. ([9]). Such 



effect is probably the highest difficulty to overc ome in both st o chastic mode l ing ( [ Sawford I 
( 200l[ )) and theory of Lagrangian turbulence (jYakhotl (|2009[) ; IZvbin et all (|2008l )). 

We acknowledge very useful discussions with M. Cencini, A. Lanotte, L. Kadanoff and 
V. Yakhot. L.B., R.B. and F.T. thank the Flash Center at the University of Chicago for 
hospitality when part of this work was prepared. We also want to thank B. Gallagher for 
helping in the data analysis and data rendering. 
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